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A new investigation of the coexistence and competition of ferroelectricity and superconductivity 
is reported. In particular we show that the starting Hamiltonian of a previous study by Birman and 
Weger (2001) can be exactly diagonalized. The result differs significantly from mean-field theory. A 
Hamiltonian with a different realization of the coupling between ferroelectricity and superconduc- 
tivity is proposed. We report the results for mean-field theory applied to this Hamiltonian. We find 
that the order parameters are strongly affected by this coupling. 
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In the present paper we report two related results from 
a reexamination of previous work on nearly ferroelectric 
superconductors 1 . In that paper a coupling term was in- 
troduced into the original Hamiltonian (Eq. 27 of Ref.QJ). 
An investigation of that coupling term shows that it only 
gives a squeezing of the phonons and no coupling to the 
electronic pairs, and therefore the Hamiltonian can be 
diagonalized exactly. However, the double mean-field 
approximation does result in an effective coupling be- 
tween the two subsystems. (Eq. 34 of Ref. Q). Therefore 
the results of the analysis of that equation remain valid. 
Our second result follows from introducing a different bi- 
quadratic coupling term in the Hamiltonian of Eq. 27 of 
Ref.Q, which satisfies gauge and inversion symmetries. 
This term does couple the two subsystems. We will treat 
this new Hamiltonian in mean-field approximation. 

We start with the model of coexistence of supercon- 
ductivity and ferroelectricity proposed by Birman and 
Weger 1 . For convenience we include its main features 
here. Birman and Weger start with two separate Hamil- 
tonians for the superconducting and ferroelectric sectors. 
The superconducting sector is a mean-field reduced BCS 
pseudo-spin model (we corrected a misprint in Ref. Q) 
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where Ek is the single-electron energy and is the pair- 
ing interaction energy. The pseudo-spin operators j pk 
obey SU (2) commutation relations and are defined as 
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where bk and bt are pair operators defined by 
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and fifef and a k ^ are the electron annihilation and creation 
operators for wave vector fc, spin ("J") , etc. 



The Hamiltonian of the ferroelectric sector of the 
model is simply a Hamiltonian of a displaced harmonic 
oscillator 
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where Bq , Bg have boson commutation relations. 

Combining those two sectors together and introducing 
a gauge invariant coupling between them, Birman and 
Weger obtain the following Hamiltonian (Eq. 27 of Ref.Q) 
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If a single fc mode is isolated 
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Here we notice that the structure of the pseudo-spin 
operators in terms of the pair operators implies 
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Jl = J2 = J3 



and therefore jf, j 2 an d J3 are invariant to any rotation 

of the form exp n ■ j 

We use this property to diagonalize the Hamiltonian 
analytically by applying the following unitary transfor- 
mation 
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= H D + j 2 (2e sin 6 - 4A cos 6) + 
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Here Ho is a Hamiltonian diagonal in the Nb) basis: 
H D = 2j 3 (2A sin 9 - e cos 9) + 

+ ujq (n b ch 2/3 + sh 2 13 + ~ + ae" 2 " 
+ 2a 7 ie-' 3 + 72i|e" 2/3 (2JV B + 1 + 4a : 

and if we want the non-diagonal parts to vanish we re- 
quire 

2e sin 9 = 4 A cos 6 

uj ae^ 2(} + 7ie~ /3 + 4a7 2 j 2 e~ 2/3 = 
72j1e~ 2 ^^sh2/3 

These yield the following set of coupled equations for the 
parameters 
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Inserting them back into the transformed Hamiltonian 
gives 
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and this is the decoupled Hamiltonian of a squeezed and 
shifted harmonic oscillator and a spin of one-half in a 
magnetic field. The energy spectrum is given by 
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We can now use the exact eigenstates of the Hamilto- 
nian to compute the order parameters of the system: 
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While the order parameter of the ferroelectric sector is 
affected by the SC-FE coupling strength, the order pa- 
rameter of the superconducting sector is not. This oc- 
curs since in the coupling term 72J 2 (j^l + Bo) , jl is a 
c-number and not an operator. 

Using a double mean-field approximation for the 
Hamiltonian (Eq. 34 of Ref. Q) and variational coherent 
state wave function, Birman and Wegeri obtained the 
following relations for the order parameters 
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with 
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Here both order parameters are affected by the SC-FE 
coupling strength (72) . 

To sum up this part, we have obtained an exact an- 
alytical solution of the ferroelectricity and superconduc- 
tivity coexistence model, used in Ref.QJ We showed here 
that contrary to the mean-field results, an exact solution 
of the model demonstrates that the ferroelectric order 
parameter is affected by the coupling, but the supercon- 
ducting order parameter is not. If one were to initially 
model the SC-FE system by Eq. 34 of Ref. [lj, the results 
reported in Ref. Q would remain valid. 

Returning to the original question of the proper cou- 
pling term in a quantum mechanical Hamiltonian, we 
need to respect gauge and inversion symmetry. Since the 
superconducting gap parameter is a complex quantity the 
added term should correspond to the bilinear |A| 2 . We 
now take this as J+J-. The coupling term for the ferro- 
electric polarization, which respects inversion symmetry 

will correspond to P 2 and can be taken as ( Bq + Bo 

Thus the bilinear coupling becomes 2+3- (Bq + Bo 

We note that = j\ + f 2 + j 3 = | + % (Using j_ j + 

would have changed J3 into — J3). The coupling becomes 

72 (j3 + \) (J3q + Bo) ■ The Hamiltonian is now 
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The interaction term in this Hamiltonian differs from the 
corresponding term (Eq. 27) in Ref. Q]. This Hamiltonian 
is not exactly solvable, therefore we make a mean-field 
approximation in the form 
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Inserting Eq. [3] into Eq. 0] results in a solvable bilinear 
Hamiltonian 
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FIG. 1: The absolute values of the superconducting order pa- 
rameter and the ferroelectric order parameter as a function 
of the coupling constant 72. The left side of the ordinate cor- 
responds to small negative values of 72. Here A = e — 71 = 

LU = 1. 



Here 
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This mean-field Hamiltonian, Eq. [fjl differs from Eq. 34 
in Ref . [l| and can be diagonalized using rotation, squeez- 
ing and displacement transformations; 

W = e a {Bl~B ) e ^(Bl 2 -B 2 ) e ieh 

These yield the following set of (self-consistent) coupled 
equations for the parameters a, (3 and 9 
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This set of equations can be reduced to a polynomial 
equation, which we solved numerically. Using these nu- 
merical results we investigate the behavior of the order 
parameters 
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Fig. shows the dependence of the two order parame- 
ters on the coupling coefficient between the ferroelectric 
and superconducting subsystems. There is a smooth evo- 
lution of each of the order parameters with the coupling 
coefficient 72. Both order parameters are non- vanishing 
at 72 = 0. We note that 72 is bounded by two critical 
values, beyond which there are no real solutions for the 
set of self-consistent equations (J7J). At the positive criti- 
cal value the superconducting gap parameter r/sc reaches 
its maximum, and the polarization order parameter t\fe 
reaches its minimum. For negative values of 72 the super- 
conducting order parameter vanishes very rapidly, while 
the ferroelectric order parameter sharply diverges, when 
72 approaches its negative critical value. 

In conclusion we note that our new model for SC-FE 
coexistence/competition does agree with the spirit of the 
Matthias conjecture that each of these cooperative effects 
tends to exclude or suppress the other—. 
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